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Abstract 

The phenomenon of quantum nucleation is studied in a nanometer-scale an- 
tiferromagnet with biaxial symmetry in the presence of a magnetic field at an 
arbitrary angle. Within the instanton approach, we calculate the dependence 
of the rate of quantum nucleation and the crossover temperature on the ori- 
entation and strength of the field for bulk solids and two-dimensional films of 
antiferromagnets, respectively. Our results show that the rate of quantum nu- 
cleation and the crossover temperature from thermal-to-quantum transitions 
depend on the orientation and strength of the field distinctly, which can be 
tested with the use of existing experimental techniques. 
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One of the most striking manifestations of the quantum character of nature is quantum 
tunneling of particles. Following the idea suggested by Caldeira and Leggett,! the tunneling 
of macroscopic object, known as Macroscopic Quantum Tunneling (MQT), has become one 
of the most fascinating phenomena in condensed matter physics. During the last decade, the 
problem of quantum tunneling of magnetization in nanometer-scale magnets has attracted a 
great deal of theoretical and experimental interest .11 The magnetic MQT includes quantum 
reversal of the magnetization (or the Neel) vector in small single-domain ferromagnets (or 
antiferromagnets), quantum nucleation of magnetic bubbles, quantum depinning of domain 
walls from defects in bulk magnets, and resonant spin tunneling in molecular magnetic 
clusters.! MQT in magnetic systems are interesting from a fundamental point of view as 
it can extend our understanding of the limits between quantum and classical physics. On 
the other hand, MQT is important to the reliability of small magnetic units in memory 
devices and the designing of quantum computers in the future. And the measurement of 
magnetic MQT quantities such as the tunneling rates could provide independent information 
about microscopic parameters such as the magnetocrystalline anisotropics and the exchange 
constants. All this makes magnetic quantum tunneling an exciting area for theoretical 
research and a challenging experimental problem. 

The problem of quantum nucleation of a stable phase from a metastable one in ferromag- 
netic films is an interesting fundamental problem which allows direct comparison between 
theory and experiment.! Consider a ferromagnetic film with its plane perpendicular to the 
easy axis determined by the magnetocrystalline anisotropy energy depending on the crystal 
symmetry. A magnetic field H is applied in a direction opposite to the initial easy direction 
of the magnetization M, which favors the reversal of the magnetization. The reversal occurs 
via the nucleation of a critical bubble, which then the nucleus does not collapse, but grows 
unrestrictedly in volume. If the temperature is sufficiently high, the nucleation of a bubble 
is a thermal overbarrier process, and the rate of thermal nucleation follows the Arrhenius 
law T T oc exp (—lI/hsT), with being the Boltzmann constant and U being the height of 
energy barrier. In the limit of T — > 0, the nucleation is purely quantum- mechanical and the 
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rate goes as Tq oc exp (— Sd/H), with S c i being the classical action or the WKB exponent 
which is independent of temperature. Because of the exponential dependence of the thermal 
rate on T, the temperature T c characterizing the crossover from quantum to thermal regime 
can be estimated as ksT c = hU/Sd- 

A few theoretical studies of the problem of quantum nucleation have been around for some 
time. Privorotskii estimated the exponent in the rate of quantum nucleation based on the 
dimensional analysis.! Chudnovsky and Gunther studied the quantum nucleation of a thin 
ferromagnetic film in a magnetic field along the opposite direction to the easy axis at zero 
temperature by applying the instanton method.! Later, Ferrera and Chudnovsky extended 
the quantum nucleation to a finite temperature.! Kim studied the effect of an arbitrarily 
directed magnetic field on the quantum nucleation of magnetization.ll The phenomenon of 
quantum nucleation was also found in nanometer-scale antiferromagnets, where the Neel 
vector is the tunneling entity. Theoretical studies on small single-domain antiferromagnets 
showed that quantum tunneling should show up at higher temperatures and higher fre- 
quencies than in single-domain ferromagnets of similar size.! This makes nanometer-scale 
antiferromagnets more interesting for experimental test. 

Up to now theoretical studies on quantum nucleation in antiferromagnets&^S were 
confined to the condition that the metastable state is caused by the magnetocrystalline 
anisotropy, which is not easily controlled in experiments. It is well-known that a magnetic 
field is a good external parameter to make the phenomenon of MQT observable. The pur- 
pose of this paper is to extend the previous consider at ions!~0 to a system with a magnetic 
field applied in an arbitrary direction between perpendicular and opposite to the initial easy 
axis (z axis). By applying the instanton method in the spin-coherent-state path-integral 
representation, we present the numerical results for the WKB exponent in quantum nu- 
cleation of a thin ferromagnetic film with the magnetic field applied in a range of angles 
7r/2 < 9h < 7r and Oh = ft respectively, where Oh is the angle between the initial easy axis 
(z axis) and the field. We also discuss the Oh dependence of the crossover temperature T c 
from purely quantum nucleation to thermally assisted processes. Our results show that the 
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distinct angular dependence, together with the dependence of the WKB exponent on the 
strength of the external magnetic field, may provide an independent experimental test for 
quantum nucleation in an antiferromagnetic film. 

For a spin tunneling problem, the rate of magnetization reversal by quantum tunneling 
can be determined by the imaginary-time transition amplitude from an initial state \i) to a 
final state |/) as 

Ufi = (/| e- HT \i)=Jv {M (r, r)} exp (S E /h) , (1) 
where S E is the Euclidean action which includes the Euclidean Lagrangian density C E as 

S E = I drd 3 rC E . (2) 



The system of interest is an antiferromagnet of about 5~10 nm in radius at a temperature 
well below its anisotropy gap. According to the two-sublattice model,! there is a strong 
exchange energy ni! ■ m 2 /x± between two sublattices, where m] and m 2 are the magneti- 
zation vectors of the two sublattices with large, fixed and unequal magnitudes, and x± is 
the transverse susceptibility. In the semiclassical regime and using the spin-coherent-state 
path- integral, one gets the Euclidean Lagrangian density for antiferromagnets (neglecting 
dissipation with the environment) as 



r a, i \\ .mi + m 2 fd(j)\ .m ( ' d<f>\ x± 

C E [B(r, r), 0(r, T )] = z— — j - j co s £ + — 



+ \a [{V6f + (V0) 2 sin 2 9] + E(9, 0), (3) 



where 7 is the gyromagnetic ratio, a is the exchange constant,0 and r = it is the imaginary- 
time variable. The E(6, <p) term includes the magnetocrystalline anisotropy and the Zeeman 
energies, m = mi — m 2 = frys, where s is the excess spin due to the noncompensation of 
two sublattices. The polar coordinate 9 and the azimuthal coordinate <fi in the spherical 
coordinate system with 1 • z = cos 9, 1 is the Neel vector of unit length and z is a unit 
vector along the z axis. The first term in Eq. (3) is a total imaginary-time derivative, 



which has no effect on the classical equations of motion, but it is crucial for the spin-parity 
effects. fffiiHl! However, for the closed instanton trajectory described in this paper (as shown 
in the following), this time derivative gives a zero contribution to the path integral, and 
therefore can be omitted. In the semiclassical limit, the rate of quantum nucleation, with 
an exponential accuracy, is given by 

T Q cx exp [Sr/H] , (4) 

where S^ m is obtained along the trajectory that minimizes the Euclidean action Se- 

In this paper, we study the quantum nucleation of the Neel vector in antiferromagnets 
with biaxial symmetry in the presence of a magnetic field at arbitrary angles in the ZX 
plane, which has the following magnetocrystalline anisotropy energy 

E (9, 0) = K\ sin 2 9 + K 2 sin 2 9 sin 2 — mH x sin 9 cos — mH z cos 9, (5) 

where K\ and K 2 are the longitudinal and the transverse anisotropy coefficients respectively, 
and K\ > 0. In the absence of the magnetic field, the easy axes of this system are ±z for 
K\ > 0. And the field is applied in the ZX plane at tt/2 < 9^ < tt. By introducing the 
dimensionless parameters as 

K 2 = K 2 /2K h H x = H x /H , H z = H z /H , (6) 

Eq. (5) can be rewritten as 

E (9, 4>) = - sin 2 9 + K 2 sin 2 9 sin 2 <p-H x sin 9cos<p- H z cos 9 + E , (7) 
2 

where E(9,<p) = 2K\E (9, 0), and Ho = 2K\/m. At finite magnetic field, the plane given 
by = is the easy plane, on which E (9, 0) reduces to 

E(9,4> = 0) = ^sin 2 9 -Hcos(9 - 9 H ) . (8) 

We denote 9q to be the initial angle and 9 C the critical angle at which the en- 
ergy barrier vanishes when the external magnetic field is close to the critical value 
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H c (9h) (to be calculated in the following). Then, the initial angle 6 satisfies 

[dE (0,(j) = 0) /d6] g = 0, the critical angle 9 C and the dimensionless critical field H c 

satisfy both [d~E(9,(f) = 0) /d0] ^ _ = _ = and [d 2 E(6,(j) = 0) /d0 2 ] ^ _ = _ = 0, which 
leads to 

- sin (20 o ) + #sin (0 O - H ) = 0, (9a) 

X - sin (20 c ) + # c sin % - H ) = 0, (9b) 

cos (29 c ) + # c cos (0 C - H ) = 0. (9c) 

After some algebra, the dimensionless critical field if c (0h) and the critical angle C are found 
to be 

H c = [(sin H ) 2/3 + |cos 9 H \ 2/S ] , (10a) 

= 1 2 '7 t f' l 'l . (i«b) 

1 + | cot 0ij-| 7 

Then the critical field of this system is if c = if c (2K 1 /m), where m = hr/s/V, and s is the 
excess spin of antiferromagnet due to the noncompensation of two sublattices. 

It is noted that for the nanometer-scale antiferromagnet at finite magnetic field, there 
exists another field H s _f m , known as the spin- flop field, which can rotate the moments of 



sublattices away from the anisotropy axis. The spin-flop field is defined as H s j = y/2HiH ex , 
with Hi = 2Kijmi being the longitudinal anisotropy field, and H ex = 2J e2 ,/m 1 being the 
exchange field between two sublattices. mi = tvyS/V is the magnetization of one sublattice, 
where S is the sublattice spin. Typical values of parameters for the antiferromagnetic 
nanoparticle are K\ ^ 10 5 erg/cm 3 and x± ^ 10~ 5 emu/G-cm 3 . The particle radius is about 
12 nm, the sublattice spin is S = 2 x 10 5 , and the excess spin is s = 10 3 . It is easy to obtain 
that the exchange energy density between sublattices is J ex (= h 2, ~f 2 S 2 /V 2 x±) ~ 1-9 x 10 10 
erg/cm 3 , H (= 2Ki/m) w 9.0 x 10 4 G, and H sJ , ss 2.8 x 10 5 G, which shows that the 
critical field is smaller than the spin-flop field. Therefore, the small e = 1 — H/H c limit can 
be performed in calculating the rate of quantum nucleation of the Neel vector, at which the 
two-sublattice configuration is still valid for antiferromagnets at finite magnetic field. 



Now we consider the limiting case that the external applied magnetic field is slightly 
lower than the critical field, i.e., e = 1 — H/H c <C 1. At this practically interest- 
ing situation, the barrier height is low and the width is narrow, and therefore the tun- 
neling rate is large. Introducing 77 = 9 C — 9 (\rj\ <C 1 in the limit of e <C 1), expanding 
[dE (9, = 0) /d9] e=6o = about 9 C , and using the relations [dE (9, = 0) /dO\ 0=e = 
and [d 2 E(9,(p = 0) /d# 2 ] - = - = 0, Eq. (9a) becomes 

sin (29 c ) (e - - 77 cos (29 c ) (2e - rj 2 ) = 0. (11) 

Simple calculations show that 77 is of the order of y/e. Thus the order of magnitude of the 
second term in Eq. (11) is smaller than that of the first term by y/e and the value of 77 is 
determined by the first term, which leads to 77 ~ y/2e/3. However, when 9h is very close to 
7r/2 or 7r, sin (29 c ) becomes almost zero, and the first term is much smaller than the second 
term in Eq. (11). Then 77 is determined by the second term when 9h — tt/2 or 7r, which 
leads to 77 ~ y/2e for 9 H ~ ir/2 and rj ~ for ~ 7r. Since the first term in Eq. (11) is 
dominant in the range of values, 9 C , which satisfies tan (29 c ) > O (y/e), rj ~ \/2e/3 is valid 
for 7r/2 + O (y/e) < 9 H < n - O (y/e) by using Eq. (10b). Therefore, 77 ~ \/2e, 0, and 
y/2e/3 for 9 H ~ 7r/2, 7r, and 7r/2 + O (y 7 ?) < 9 H < n — O (y/e), respectively. In this case the 
potential energy E (9, 0) reduces to the following equation in the limit of small e, 

E (5, 0) = K 2 sin 2 sin 2 (9 + 5) + H x sin (9 + 5) (1 - cos 0) + ^ (5) , (12) 

where <5 = # — 9 (\5\ <^ 1 in the limit of e 1), and E\ (8) is a function of only 5 given by 

^ (5) = sin (20 c ) (35 2 ri - 5 3 ) + X - cos (29 c ) 

It can be shown that in the region of n/2 + O (y/e) < 9h < tt — O (y/e), O (y/e) < 9 C < 
7r/2 — O (y/e), rj and 5 are of the order of y/e, the second term in Eq. (13) is smaller than 
the first term in the small e limit. It is convenient to use the dimensionless variables 

r' = e 1/4 r/r , r' = e 1/4 u r, 5 = 8/yfe, (14a) 

where r = yj a/2K 1: and 
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6 ' 2 f 



+ 5 3 V - -5 A 
4 



(13) 



sin 9 C 



+ 



-1/2 



m \H x + 2K 2 sm9 c 
Then the Euclidean action Eq. (2) for 7r/2 + O (y/e) < 9 H < n — O (y/e) becomes 



(14b) 



S E p(r',r') ,0(r',r')] = ^ / dr'd'r \ ^r-wj, y 



eujQ 
4 



&r 2 



^eV^ sin2 (, + ^ (|Q _ ^ sm {9o + v^) 
Z 2 sin 2 sin 2 (0„ 4- v 7 ^) + 2if x sin 2 ^0 sin (0 O 4- v 7 ^) 



2 7 
4-2^ 



4-^ 3/2 (V5) 2 + ie 172 sin 2 (0 O 4- (V'0) 2 



4--e 



3/2 



(15) 



where A = sin(2# c ). In Eq. (15) we have performed the integration by part for the term 
— i^cos^(^) and have neglected the total imaginary-time derivative. In can be showed 
that for tc/2 + O (y/e) < 9 H < n — O (y/t), only small values of contribute to the path 
integral, so that one can replace sin 2 in Eq. (15) by 2 and neglect the term including 
(V'0) 2 which is of the order e 2 while the other terms are of the order e 3//2 . Then the Gaussian 
integration over leads to 



J V{5(r',r')}e W ("J^") 



(16) 



where the effective action is 

Sf[5(r',r')] =h 8 e 1 ' 2 r* 



sm i 



H x + 2K 2 sin9 
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+ 



2X^!^ 1/2 



x / dr'd 3 r' 



(17) 



and s = m/tv) is the excess spin of antiferromagnets. Introducing the variables r = r'y/A 
and r = r'y/A, the effective action Eq. (17) is simplified as 



S e J f \5(r,r)} = hse 1 ^ 



1 ( sm9 r 2\./m^ 2 



A \H x + 2K 2 sm6 t 

2 



+ 



(18) 



For the quantum reversal of the Neel vector in a small particle of volume V <C Tq, the 
Neel vector is uniform within the particle and S does not depend on the space f. In this 
case Eq. (18) reduces to 
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H x + 2K 2 sin 9 C m 
The corresponding classical trajectory satisfies the equation of motion 



1/2 



x I dr 



(19) 



d 2 5 1 ^ 3- 



dr 2 

Eq. (20) has the instanton solution 



^ = -V6J--7. (20) 



- _ 

5(T) = cosh 2 (3V4 X 2-V4^)' (21) 



corresponding to the variation of 5 from <5 = at r = — oo, to S = y6e at r = 0, and then 
back to 5 = at r = oo. The associated classical action is found to be 

s i = 2 i^ x3 i/4 ^ 5/4 1^ gg | V6 / l + |cotg g | 2/3 + 2 X ±^A 

5 36 y/l + |cot^| 2/3 V 1 -e + 2K 2 (l + |cot^| 2/3 ) m2 J 

(22) 

In the WKB approximation, the tunneling rate T of a particle escaping from a metastable 
state has the relation T cc exp (— B/ti). The WKB exponent B is approximately given 
by U/uJb, where U is the height of barrier, and oo^ is the frequency of small oscillations 
around the minimum of the inverted potential and characterizes the width of the barrier 
hindering the decay process. For magnetic quantum tunneling, u 2 (= -E'[ (8 m ) /Af) is 
inversely proportional to the effective mass of the magnets, where the mass is induced by the 
transverse component of magnetic field, and 5 m corresponds to the position of the minimum 
of the inverted potential. For general case, the WKB exponent B should be proportional 
to the power of the parameter e = 1 — H/H c since the height and the width of barrier are 
proportional to the power of e. Simple analysis of Eq. (8) shows that the value of e should be 



small in order to obtain a large tunneling rate. For this case, we obtain the height of barrier 
as E 1 (= U/2K 1 V) = sin (20 c ) (6e) 3/2 /27 at <5 m = 2^/3 and the oscillation frequency 
around the minimum of the inverted potential —E% (5) as Ub = 2 -1 / 4 x 3 1 / 4 e 1 / 4 u;o\/sin (29 c ), 
where ujq is shown in Eq. (14b). Then we approximately obtain the WKB exponent B as 

„ U 
B 

5 /4 |cotg g | 1/6 I l + |cotfl g | 2/3 , 2X±K\ * 

" 37/4 ^i + m^V 1 -^^^!^^) m2 ' ' 

which agrees up to the numerical factor with the result in Eq. (22) obtained by using the 
explicit instanton solution. 

For a large non-compensation (m ^> y/xU^i) , Eq. (23) reduces to 

5 ™ = fe^g (24) 

l-e + 2K 2 (l + |cot0#| 2/3 ) 

which agrees well with the result of quantum tunneling of magnetization in single-domain 
ferromagnetic particles with biaxial symmetry in an arbitrarily directed field.0 For a small 
non- compensation (m <C y/x±Ki) > Eq. (23) reduces to the result for single-domain antifer- 
romagnetic particles, 

S AFM = 2lV4 X 3V4 VxTKl y c 5/4 |cot^| 1/6 



5 " ' • ''l • COt^^ 



(25) 

which is in good agreement with the result in Refs. 20 and 21. 

Now we turn to the nonuniform problem. In case of a thin film of thickness h less than 
the size r /e 1//4 of the critical nucleus and its plane is perpendicular to the initial easy axis, 
we obtain the effective action after performing the integration over the z variable in Eq. 
(19), 



x / drcPv 



(26) 



At zero temperature the classical solution of the effective action Eq. (26) has O (3) symmetry 
in two spatial plus one imaginary time dimensions. Therefore, the solution 5 is a function of 
u, where u = (p 2 + T 2 ) 1 ^ 2 , and p = (x 2 + y 2 ) 1 ^ 2 is the normalized distance from the z axis. 
Now the effective action Eq. (26) becomes 



1/2 



(27) 



The corresponding classical trajectory satisfies the equation of motion 



^ + £^ = ^J-V. ( 28) 
du 2 udu 2 4 y ' 

By applying the similar method,&0 the instanton solution of Eq. (28) can be found numeri- 
cally and is illustrated in Fig. 1. The maximal rotation of the Neel vector is <5 max ~ 6.8499 
at T = and p = 0. Numerical integration in Eq. (27), using this solution, gives the rate of 
quantum nucleation for a thin antiferromagnetic film as 

T Q oc exp {-S E /h) 



exp { -74.39se 3/4 r^ 



2/3 



I cot e 



1/6 



1 + 


cot 9h 


2/3 


1 - e + 2K 2 




cot 6h 


2/3\ 



1/2 



+ . (29) 



m 2 



For a large non-compensation, Eq. (29) reduces to the result for quantum nucleation in a 
thin ferromagnetic film 



Y Q a exp {S F E M /ti) 



1 + I cot H 



1 2/3 



exp <^ -74.39se 3/4 r 2 /i = 1 H 1 ) . (30) 



Icot ^wl 176 J 1 - e + 2A^ 2 [ 1 + Icot ^ ff | 2/3 



For a small non- compensation, Eq. (29) reduces to the result for quantum nucleation in a 
thin antiferromagnetic film 
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T Q a exp {-Sf M /h) 



= exp J -105.2^e3/V ^V 1 + |COt ;f /3 I . (31) 
( 7 |cot^| 1/6 J 

At high temperature, the nucleation of the Neel vector is due to thermal activation, 
and the rate of nucleation follows IV oc exp (— W m i n /kBT), where W m i n is the minimal 
work necessary to produce a nucleus capable of growing. In this case the instanton solution 
becomes independent of the imaginary-time variable f. In order to obtain W min , we consider 
the effective potential of the system 

U eff = J dh [| ((W) 2 + sin 2 9 (V0) 2 ) + E (9, 0)] . (32) 

For a cylindrical bubble Eq. (32) becomes 



From the saddle point of the functional the shape of the critical nucleus satisfies 



(33) 



dp 2 + pdp~ 2 I' 



3=2 + =3= = V*-7*- (34) 



The solution can be found by numerical method similar to the one in Refs. 5 and 7. Fig. 2 
shows the shape of the critical bubble in thermal nucleation, and the maximal size is 3.906 
at p — 0. Using this result, the minimal work corresponding the thermal nucleation is 

Wmin = 41.3376^ier 2 /i. (35) 

Comparing this with Eq. (29), we obtain the approximate formula for the temperature 
characterizing the crossover from thermal to quantum nucleation as 

-1/2 



, . cot g H 1/6 / 1+ cot 9 H 2/3 2x±iTi , 

k B T c « 0.55— — , — Y + -—5- | . (36) 

/ l + |cot^| 2/3 \l-e + 2K 2 [l + cot^ 2/3 J 



m 2 



For the FM case, i.e., the case of large non-compensation, the crossover temperature is 



:F.w , „ rr ^ie 1/4 |cot^| 1/6 



k B T b c M « 0.55— . 1 (37) 



1-6+2^9 1+ Icotfl 



7 # 



1 2/3 
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While for the AFM case, i.e., the case of small non-compensation, 



k B T. 



^ M «0.39fryJ— e 1/4 

V x± 



I cot o H 1 1/6 



(38) 



1 + |cot^| 2/3 

To observe the quantum nucleation one needs a large crossover temperature and not too 
small a nucleation rate. Note that x± = m 1/Jex and mi = hr/S where J ex is the exchange 
interaction between two sublattices and S is the sublattice spin, Eq. (36) can be written as 

.Kxe^hrf | cot H | 1/6 



k B T c « 0.55- 



m 



I cot 0# 



2/3 



X 



1 + 


cot Oh 


2/3 


1 - e + 2Z 2 


(i + i 


cot 


2/3\ 




-1/2 



In Fig. 3, we plot the 0# dependence of the crossover temperature T c for typical values of 
parameters for nanometer-scale antiferromagnets: K\ = 10 7 erg/cm 3 , J ex = 10 10 erg/cm 3 , 
K 2 = 1, m = 10 emu/cm 3 , S/s — 100, e = 0.01 in a wide range of angles 7r/2 < 9h < 7T. 
Fig. 3 shows that the maximal value of T c is about 1.916K at = 2.350, which is one 
or two orders of magnitude higher than that for ferromagnets with a similar size.iQ Note 
that, even for e as small as 10 -3 , the angle corresponding to an appreciable change of the 
orientation of the Neel vector by quantum tunneling is <5 2 = \/6e rad> 4°. The maximal 
value of T c as well as Tq is expected to be observed in experiment. 

Now we study the situation that the magnetic field is applied opposite to the initial easy 
axis, i.e., 6h = 7r. In this case, 0q = C = 0, H x = 0, and rj = from Eqs. (10a) and 
(10b). By using the dimensional variables r' = e 1 / 2 r/r , r' = ^^ujqT, 5 = where 
r = y/a/2K±, and 



UJ 



m 



1 | 2 X ±K 1 



2K 2 



we obtain the Euclidean action as 
^p(r',r'), </>(r',r')] 



e 2 aj 



'J3 / J X-L 5, ,2 



2r 



e u. 



(^) +f [d? 



12 



■ m 3/2 

7 



-f2 



(39) 



The Gaussian integration over reduces Eq. (39) to the following effective action 

— 2 - - ' r 
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sy> p ( ,,, )H 2M/«v [I ( * )* + i (v^ + i - ? 



(40) 



In the case of quantum nucleation of the Neel vector in a small particle of volume V <C Tq, 
the Neel vector is uniform within the particle and Eq. (40) reduces to 



S e J f [5(r')} 



2K X 



e 3 /V / dr 



lfd5\ 2 I f-2 5^ 
2{d?) + 2{ 5 "4 



(41) 



The classical trajectory satisfies the equation of motion 

d 2 S , l-=a 



which has the instanton solution 



5(r') 



(42) 



(43) 



cosh t' ' 

corresponding to the variation of 5 from 5 = at r = — oo, to 5 = 2^/e at r = 0, and then 
back to 5 = at r = oo. Substituting this solution into Eq. (41) we obtain that 



5d = 5 fisf3/v M + 2x±M 

3 \2K 2 m 2 J 



1/2 



(44) 



In the case of quantum nucleation in a thin film of thickness h less than the size r /y/e of 
the critical nucleus with its plane perpendicular to the initial easy axis, the Euclidean action 
Eq. (39) becomes 



S e J f [5 (r', r')] = Anhe^rlh ( JL + 2 -^^j ^ J dun 



1 fd5\ 2 1 (-2 5^ 
2\Tu) + 2{ 6 ~ T 



(45) 



where u 2 = p' 2 + r' 2 , and p' 2 = x' 2 + y' 2 . The corresponding classical equation of motion 
satisfies 
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The instanton solution of Eq. (46) can be found numerically and is illustrated in Fig. 4. 
Numerical integration in Eq. (45), using this solution, gives the WKB exponent for the 
subbarrier bubble nucleation in an antiferromagnetic film as 



T Q oc exp (S E /H) = exp | -37.797se 1/2 r^ ^ 



-+- 

2 



2K 2 m 2 



1/2' 



(47) 



For a cylindrical bubble, the effective potential of the thermal nucleation is found to be 



POO 

U eff = AnKierlh I dp' p 1 
Jo 



1 fd5\ 2 1 [-2 5^ 
2\W + 2V "T 



(48) 



wherein the shape of the critical nucleus corresponds to a saddle point of this functional: 

d 2 5 1 dS — 1-3 , An . 

Eq. (49) can be solved by the numerical approach, and the solution is showed in Fig. 
5. Numerical integration of this solution in Eq. (48) gives T T oc exp (— W min /k B T) with 
W^min = 23.402^16^/1. Then the crossover temperature is found to be 

K x e x l A ( 1 2 X ±K i y 1/2 



kBT ^ 0M9 ^_ 1^;+^] (50) 



In conclusion, we have investigated the quantum nucleation of the Neel vector in 
nanometer-scale antiferromagnets with biaxial symmetry in the presence of an external mag- 
netic field at arbitrary angle. By applying the instanton method in the spin-coherent-state 
path-integral representation, we obtain the analytical formulas for quantum reversal of the 
Neel vector in small magnets and the numerical formulas for quantum nucleation in thin 
antiferromagnetic film in a wide range of angles n/2 < 6h < 7r, and 6h = ^ respectively. 
The temperature characterizing the crossover from the quantum to thermal nucleation is 
clearly shown for each case. Our results show that the rate of quantum nucleation and 
the crossover temperature depend on the orientation of the external magnetic field dis- 
tinctly. When 6 H = ir, the magnetic field is applied antiparallel to the anisotropy axis. It 
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is found that even a very small misalignment of the field with the above orientation can 
completely change the results of tunneling rates. Another interesting conclusion concerns 
the field strength dependence of the WKB exponent or the classical action in the rate of 
quantum nucleation. We have found that in a wide range of angles, the e (= 1 — H/H c ^j 
dependence of the WKB exponent or the classical action S c i is given by e 5//4 , not e 3 / 2 for 
9jj = 7i. Therefore, both the orientation and the strength of the external magnetic field 
are the controllable parameters for the experimental test of quantum nucleation of the Neel 
vector in nanometer-scale antiferromagnets. If the experiment is to be performed, there are 
three control parameters for comparison with theory: the angle of the external magnetic 
field On-, the strength of the field in terms of e, and the temperature T. 

Recently, Wernsdorfer and co-workers have performed the switching field measurements 
on individual ferrimagnetic and insulating BaFeCoTiO nanoparticles containing about 10 5 - 
10 6 spins at very low temperatures (0.1-6K).@ They found that above 0.4K, the magneti- 
zation reversal of these particles is unambiguously described by the Neel-Brown theory of 
thermal activated rotation of the particle's moment over a well defined anisotropy energy 
barrier. Below 0.4K, strong deviations from this model are evidenced which are quanti- 
tatively in agreement with the predictions of the MQT theory without dissipation. Em it 
is noted that the observation of quantum nucleation of ferromagnetic or antiferromagnetic 
bubbles would be extremely interesting as the next example, after single-domain nanopar- 
ticles, of macroscopic quantum tunneling. The experimental procedures on single-domain 
ferromagnetic nanoparticles of Barium ferrite with uniaxial symmetry!! may be applied to 
the antiferromagnetic systems. Note that the inverse of the WKB exponent B" 1 is the 
magnetic viscosity S at the quantum-tunneling-dominated regime T C T c studied by mag- 
netic relaxation measurements.^ Therefore, the quantum nucleation of the antiferromagnetic 
bubbles should be checked at any 9h by magnetic relaxation measurements. Over the past 
years a lot of experimental and theoretical works were performed on the spin tunneling in 
molecular Mn 12 - and Fe 8 @ clusters having a collective spin state S = 10 (in this paper 
S = 10 6 ). Recent progresses in the field of molecular magnets include the measurements of 
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the spin-lattice relaxiation rate,El the specific heat and the ac-susceptibility on Mni2-Ac,l 



the studies of energy sublevels of the ground state by inelastic neutron scattering 



and the zero-field magnetic relaxation of Mn^-AcJ 3 ^ and the model calculation of the mag- 
netization relaxation based on phonon-assisted tunneling on Mni2-Ac;0 the INS experiment 
on Fe§,ii the measurement of effects of nuclear spins on the magnetic relaxation of Feg,S 
and the nonadiabatic Landau-Zener tunneling in Fesl^S Further experiments should focus 
on the level quantization of collective spin states of S = 10 2 -10 4 . 

The ferromagnet or antiferromagnet is typically an insulating particle with as many as 
10 3 ~ 10 6 magnetic moments. For the dynamical process, it is important to include the effect 
of the environment on quantum tunneling caused by phonons,@0 nucleation spins,^! and 
Stoner excitations and eddy currents in metallic magnets.il However, many studies showed 
that even though these couplings might be crucial in macroscopic quantum coherence, they 
are not strong enough to make the quantum tunneling unobservable.Bil~i! The theoretical 
calculations performed in this paper can be extended to the antiferromagnetic systems with 
more general symmetries, such as trigonal, tetragonal and hexagonal symmetries. We hope 
that the theoretical results presented in this paper may stimulate more experiments whose 
aim is observing quantum nucleation in nanometer-scale magnets. 
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Figure Captions: 

Fig. 1 The instanton, corresponding to subbrrier bubble formation in a thin film by quan- 
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turn tunneling in a magnetic field with n/2 < Oh < ft, for r = 0, r = ±0.5, r = ±1, 
T = ±1.5, and r = ±2. 

Fig. 2 The shape of the critical bubble in a thermal nucleation of magnetization in a 
magnetic field with ir/2 < 9 H < ft- 

Fig. 3 The 0h dependence of the crossover temperature T c for ir/2 < 6h < ft- Here, 
Ki = 10 7 erg/cm 3 , J ex = 10 10 erg/cm 3 , K 2 — 1, m — 10 emu/cm 3 , S/s — 100, and 
e = 0.01. 

Fig. 4 The instanton, corresponding to subbrrier bubble formation in a thin film by quan- 
tum tunneling in a magnetic field with 6h — ft, for r = 0, r = ±0.5, f = ±1, r = ±1.5, 
and r = ±2. 

Fig. 5 The shape of the critical bubble in a thermal nucleation of magnetization in a 
magnetic field with Oh = ft- 
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